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Abstract 

<N ■ We investigate non-Gaussian isocurvature perturbations generated by the evo- 

lution of Goldstone modes during inflation. If a global symmetry is broken before 



inflation, the resulting Goldstone modes are disordered during inflation in a precise 



and predictable way. After inflation these Goldstone modes order themselves in a self- 
similar way, much as Goldstone modes in field ordering scenarios based on the Kibble 



mechanism. For (Hf n j / M^) ~ 10 6 , through their gravitational interaction these 



Goldstone modes generate density perturbations of approximately the right magni- 
tude to explain the cosmic microwave background (CMB) anisotropy and seed the 
structure seen in the universe today. We point out that for the pattern of symmetry 
breaking in which a global U(l) is completely broken, the inflationary evolution of 
the Goldstone field may be treated as that of a massless scalar field. Unlike the more 
commonly discussed case in which a global U(l) is completely broken in a cosmolog- 
ical phase transition, in the inflationary case the production of defects can be made 
^ ■ exponentially small, so that Goldstone field evolution is completely linear. In such a 

model non-Gaussian perturbations result because to lowest order density perturba- 
tions are sourced by products of Gaussian fields. Consequently, in this non-Gaussian 
model N-point correlations may be calculated by evaluating Feynman diagrams. We 
explore the issue of phase dispersion and conclude that this non-Gaussian model 
predicts Doppler peaks in the CMB anisotropy. 
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1. Introduction 



A key problem of modern cosmology is the origin of the primordial fluctuations 
that later led to structure in the universe. From a theoretical perspective, at this 
point the most satisfying theory for the very early universe is inflationary cosmol- 
ogy because of its nice resolution of the horizon, smoothness, flatness, and monopole 
problems! 11 In its simplest form, inflation with a single real scalar field predicts an 
approximately scale-invariant spectrum of Gaussian adiabatic primordial density per- 
turbations! 21 But if other light scalar fields are present during inflation (light relative 
to the expansion rate i?j n / during inflation), other types of density perturbations 
are generated as well, and these other types of perturbations leave a different and 
distinguishable imprint on the universe today. 

The observation that inflation can produce perturbations other than the com- 
monly discussed almost scale-free Gaussian adiabatic density perturbations is not 
new. Possibilities have been pointed out by various authors. In inflationary models 
with a Peccei-Quinn symmetry broken before inflation, the inflationary expansion 
disorders the axion field, which during inflation may be regarded as a massless scalar 
field. After inflation, when the axion mass becomes relevant, fluctuations in the ax- 
ion field translate into perturbations in the axion mass density. Allen, Grinstein, 
and Wise observed that a fourth-order derivative term in the effective theory for 
the Peccei-Quinn symmetry breaking gives rise to non-Gaussian fluctuations in the 
axion density. There are other possibilities for non-Gaussianity in axion density per- 
turbations from inflation; W however, in many cases inflationary axion perturbations 
are very nearly Gaussian. It has also been pointed out that for inflation with mul- 
tiple scalar fields possibilities exist for non-Gaussian fluctuations! ' 61 For example, 
when the classical slow-roll path bifurcates into two paths toward the minimum, in 
what may be described as hitting a mogul, quantum fluctuations about the slow-roll 
path determine which side of the mogul the inflaton field chooses and non-Gaussian 
fluctuations result. Fluctuations from finite size bubbles in extended inflation give 
non-Gaussian density perturbations. Also by combining inflation with various sorts 
of topological defect models non-Gaussian perturbations may be produced.' 8 ' 91 

In this paper we investigate the evolution of Goldstone modes during inflation 
and the density perturbations subsequently generated. We assume an exactly massless 
Goldstone mode and also that the potential for the inflaton field is such that the usual 



1 



adiabatic Gaussian perturbations from quantum fluctuations of the inflaton field are 
so small as to be observationally irrelevant. We assume that unlike in axion models the 
global symmetry giving rise to the Goldstone mode remains exact today. The effect 
of even a small Goldstone mass drastically changes the analysis of the perturbations. 
Vilenkin in 1982 pointed out how the ordering of an initially disordered Goldstone 
field can generate density perturbations of the right order to explain the origin of 
structure in the universe.' 101 Vilenkin considered a Goldstone field initially disordered 
through a symmetry-breaking phase transition, so that on large scale there are no 
correlations. The evolution of Goldstone modes for large N have been studied by 
Turok and Spergel.' 111 Here we assume that the global symmetry is broken before 
inflation so that the inflationary dynamics of the Goldstone fields determine the state 
of the Goldstone field on superhorizon scales at the end of inflation. 

Miiller and Schmid have also considered density fluctuations arising from an infla- 

[12] 

tionary Goldstone mode. Their analysis differs from ours in that they compute only 
correlations of the Goldstone field contribution to the stress-energy. The cosmological 
perturbations that we observe today are perturbations in the gravitational potential 
and in the total stress-energy. Since the contribution of the Goldstone modes to 
the total stress-energy is small, it is justified to ignore the effect of gravity on the 
evolution of the Goldstone modes, treating them as a stiff source. The rest of the 
matter in the universe, which gives the overwhelming contribution to the total stress- 
energy, may be treated as a single-component perfect fluid, first with the equation of 
state of radiation and later with a pressureless equation of state. The Goldstone stiff 
source excites perturbations of this single-component fluid coupled to gravity. Cos- 
mological perturbations observed today are computed using Green functions — their 
correlations are related by an integral transform to correlations at earlier times of the 
sort considered by Miiller and Schmid. 

Inflationary Goldstone modes are interesting for two reasons. First, if we are 
to confront inflation with observation, we must fully explore the possibilities offered 
by inflationary cosmology rather than exclusively focusing on the predictions of the 
simplest models. The second motivation has to do with the simplicity of this model 
and its qualitative similarity to field ordering models based on the Kibble mechanism 
(i.e., cosmic strings, global monopoles, textures, nonlinear Goldstone modes, etc.). 
For inflationary Goldstone modes, as we shall show, in many situations their evolution 
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is linear, or almost linear. The resulting density fluctuations, however, are non- 
Gaussian because they are sourced by the stress-energy from the Goldstone modes, 
which is the sum of products of two Gaussian fields rather than a Gaussian field itself. 
Because of this, the expectation values of the CMB multipole moments are four-point 
functions in a Gaussian theory, and thus comparatively simple to evaluate. Despite 
the underlying Gaussian character of these linear inflationary modes, this model is 
remarkably similar to field ordering models based on the Kibble mechanism. In both 
cases one starts (either after inflation or after the symmetry-breaking phase transition) 
with a disordered order parameter field that orders itself in a self-similar way, with the 
characteristic length scale at any time equal to the Hubble length at that time. In both 
cases the matter density perturbations are sourced by a Goldstone mode stress-energy 
0^(x, t) that is non-Gaussian, so that a simple decomposition into momentum modes 
is not possible. The principal difference between inflationary Goldstone modes and 
Goldstone modes from the Kibble mechanism lies in the correlations of the initial 
field configurations on superhorizon scales. With the Kibble mechanism there are 
initially no correlations on large scales because of considerations of causality, but in 
the inflationary case at reheating there exist correlations on arbitrarily large scales. 
However, this difference is expected to have a small effect, probably only changing 
the parameters of a scaling solution but not its qualitative behavior. 

Recently there has been much discussion in the literature about the structure of 
the small-angle CMB multipole moments in field ordering theories.' 13 171 An impor- 
tant qualitative issue is whether field ordering theories predict oscillations, or Doppler 
peaks, in the large-£ CMB multipole moments, similar to those predicted for adia- 
batic Gaussian fluctuations from inflation. For adiabatic Gaussian fluctuations each 
mode in momentum space has a definite fixed phase, although its real amplitude is a 
Gaussian random variable. Therefore modes of wavelength comparable to or smaller 
than the Hubble length at last scattering undergo acoustic oscillations with a definite, 
fixed phase. By contrast, for non-Gaussian isocurvature models these acoustic modes 
have at least some dispersion in their phase. Since these modes are sourced primarily 
during horizon crossing over a time of order the Hubble time, there is no reason to 
expect what would be the decaying mode at earlier times to be absent, and since 
the source is nonlinear, there is no reason to expect that the time dependence of the 
source for a given mode, as opposed to its overall amplitude, not to fluctuate. A fixed 
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time dependence would eliminate any phase dispersion. Therefore, it is clear that 
oscillations of the CMB moments in field ordering theories are at least to some extent 
suppressed, but how much is a detailed quantitative question. While there definitely 
exist effects that tend to smear the phase, it is not a priori clear whether these ef- 
fects are strong enough to suppress the oscillatory character of the small-angle CMB 
multipole moments. Crittenden and Turok' 141 have found oscillatory Doppler peaks 
for global textures (where the order parameter space is SU(2) = S 3 ). Their model 
for computing the CMB moments assumes a fixed time dependence for the source 

[15] 

for each mode, thus not allowing for any phase dispersion. Durrer et al. have also 
studied the CMB moments for textures, obtaining results essentially in agreement 
with Crittenden and Turok. Albrecht et al., [16] on the other hand, have suggested that 
phase dispersion may be very significant, especially for cosmic string models, which 
they claim exhibit no secondary Doppler peaks. It is hoped that the CMB anisotropy 
moments from linearized inflationary Goldstone modes, which because of its simplic- 
ity can be solved exactly in terms of four-point functions in a Gaussian theory, may 
shed some light on these questions. 

Goldstone modes arise through the breaking of a global symmetry G to a smaller 
group H, giving rise to massless modes equal in number to the difference in dimen- 
sionality of the two groups. Most studied has been the evolution of such modes after 
a cosmological phase transition. In these field ordering or topological defect models, 
it is typically assumed that initially (either because of preferred initial conditions or 
because of a prior epoch of inflation) one starts with an exactly homogeneous and 
isotropic universe at a high temperature in which the symmetry G is unbroken. Later, 
as the universe cools, a phase transition takes place in which the symmetry G breaks 
to H. Immediately after the phase transition at each point an orientation of G/H is 
chosen at random, but beyond some correlation length, which by considerations of 
causality cannot exceed the horizon size, the orientations are uncorrelated. Subse- 
quently the order parameter field orders itself in a self-similar way, described by a 
scaling solution, which in general must be determined numerically. 

In this paper we consider a situation in which the symmetry G is broken to 
H before inflation. In this case the state of the field G/H at the end of inflation is 
completely determined by the inflationary dynamics of the field G/H during inflation. 
All vestiges of initial conditions before inflation are erased. After inflation the coset 
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field evolves classically, becoming aligned on increasingly larger scales, just as in the 
scenarios based on the Kibble mechanism. Although the field ordering dynamics 
in the two cases (after the end of inflation or after the phase transition) are the 
same, the initially conditions at the respective times are quite different. After the 
phase transition the coset field is uncorrelated over large distances — one starts with a 
white noise spectrum on large scales. By contrast, at the end of inflation there exist 
correlations on arbitrarily large scales. 

In the case where Hi n f <C f g (where Hi n f is the Hubble constant during inflation 
and f g is the symmetry breaking scale, roughly the radius of the coset space), one may 
adopt the following stochastic picture of the inflationary dynamics of the coset field, 

[18] 

taken from the stochastic approach to chaotic inflation. We shall ignore fluctuations 
of the coset field on scales less than the Hubble length H~ n j, treating the coset field as 

classical and constant over a Hubble volume, a cube of volume H^. As the universe 
expands by a factor of two, each such cube subdivides into eight subcubes, each of the 
original size, and in each of these subcubes the coset field takes random step of order 
A<fi ~ H in f. The random steps in different subcubes are uncorrelated. This process 
repeats until the end of inflation. Physically, the random steps represent the freezing 
in of quantum fluctuations on subhorizon scales by the inflationary expansion. 

In the most general case, owing to nonlinearity, the evolution of the nonlin- 
ear sigma model G/H during and after inflation is quite complicated and thus not 
amenable to an analytic treatment. However in special cases nonlinearity plays a 
negligible role. When G/H = S 1 , the nontrivial topology of S 1 is irrelevant, and 
we may replace S 1 by the real line, treating the Goldstone mode as a free massless 
real scalar field. If the potential in the radial direction is sufficiently stiff, the process 
by which the field jumps over the origin creating a configuration of nonzero winding 
number is exponentially suppressed. This process may be described as the nucle- 
ation of a cosmic string loop. For the loop not to recollapse, it must have a radius 
larger than a Hubble length. Note that for this pattern of symmetry breaking, in 
scenarios based on the Kibble mechanism nonlinearity is always important, because 
there is no way to suppress the formation of cosmic strings. The other case in which 
inflationary Goldstone modes may be linearized arises when Hi n f <C f g - In this case, 
provided one does not consider too large a volume, the Goldstone field never wanders 
far enough from its average value for nonlinear effects to become important. 
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The organization of this paper is as follows. In section II we present the initial 
condition for a linearized Goldstone mode at the end of inflation and give its evolution 
into the radiation-dominated era as well as the resulting contribution of the Goldstone 
mode to the stress-energy. In section III the resulting stress-energy is renormalized 
and the effect of this source through gravity on the radiation fluid is computed. A 
derivation of the Green's functions for computing the perturbations of the radiation 
fluid as an integral transform over the stiff source stress-energy is relegated to an 
appendix. Because our computation does not include the decaying modes of the 
Goldstone field, our calculations become unreliable on subhorizon scales. We follow 
the evolution of the radiation fluid modes until horizon crossing and compute the 
phase dispersion in terms of a covariance matrix of two point functions of the two 
radiation fluid modes. We find the amount of phase dispersion to be small. In Section 
IV we present some concluding remarks indicating future directions. 



For simplicity we consider a single Goldstone mode represented by a real mass- 
less scalar field, as in the G/H = S 1 model. The generalization to several Goldstone 
modes is straightforward. After the end of inflation we may treat this field as a Gaus- 
sian ensemble of classical fields that evolve classically, just as one typically assumes 
in the discussion of the usual adiabatic perturbations from inflation. Initially, on 
superhorizon scales one has 



sian random variables subject to the constraint a(k) = a(— k)* with the two-point 
function 



2. Goldstone Mode Evolution and Stress-Energy 




(2.1) 



where i?j n / is the Hubble constant during inflation and the a(k)'s are complex Gaus- 




(2.2) 



The subsequent evolution of the Goldstone modes is given by 




(2.3) 
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where T(r), k) satisfies 



f + ^lf + k 2 T = (2.4) 



and the boundary conditions T(rj — 0, k) — 1, T(rj — 0, k) — 0. Here the spacetime 
metric is ds 2 = a 2 (if) ■ [—drj 2 + dx 2 ] and the dots denote derivatives with respect to 
conformal time rj. 

The stress-energy from the massless scalar field is given by 

= [d u <t>){dv<t>) - \guA9 aP {d*<t>){dp<t>)\- (2.5) 

Assuming a radiation-dominated universe so that a(r)) = rj, we may solve eqn. 
(2.4) explicity, obtaining 

T( m k)=j (kr)) = -^-. (2.6) 

Here we have given an explicit analytical solution for the evolution of the stiff 
source. In most field ordering theories (such as textures, global monopoles, and 
cosmic strings) the evolution of the stiff source must be computed numerically on a 
huge lattice. 

In computing matrix elements involving G M ^(k) it is necessary to renormalize. In 
this paper we are concerned with computing two-point functions of the form 

(e^(k)0 MV (-k)>. (2.7) 

These can be renormalized a la Zimmerman by replacing the unrenormalized inte- 
grand Iunrenormfe) with J renorm (k) = Iunrenormfe) lunrenormfe = 0). Given the 

approximations used here, in which the decaying mode is ignored resulting in a de- 
scription for the behavior of the scalar field accurate only on superhorizon scales, 
the need to renormalize is not completely manifest, because integrating over the un- 
renormalized integrands gives convergent integrals, with dominant contribution from 
momenta of order k m H{ n j. But k Hi n f is precisely where the approximations 
used here break down. The ultraviolet divergence has been avoided because we have 
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extrapolated the infrared enhanced correlations on superhorizon scales to subhori- 
zon scales. In a sense through our approximations we have introduced a cutoff at 
^internal ~ Hinj • The renormalized matrix elements can be computed straightfor- 
wardly by subtracting the integrand at k — 0. After this subtraction at any given 
time 77 the dominant modes contributing to are peaked around k ~ r]" 1 (i.e., of 
wavelength of order the horizon size, or Hubble length), as one might expect. 
For future reference we give the expansion of the Fourier components 



0^(k) = J d 3 x e"* k ' x 0^(x) 



(2.8) 



in terms of the Gaussian random variables a(k) and decompose these Fourier com- 
ponents into their scalar, vector, and tensor parts. 

Substituting eqn. (2.3) into eqn. (2.5), one obtains 



0oo(k) = ^/ /' (H < «(k-k>(k') 



x 



(2tt) 3 |k-k'|§ |k'|§ 
l dT( V ,\k-k'\)dT( V ,\k'\) 1 

2 aj ft) 2 (k " k) ' k n^lk-kDT^Ikl) 



(2.9) 



©o,(k) = ^ | 



d s k' a(k-k>(k') 
|k-k'|§ |k'|§ 



9ry 



Tfa, |k'|) + i (k - k^ g^Jgl) T (^, |k - k'| 



(2.10) 



e y (k) = 



r d 3 k' a(k-k'Wk') 



inf 



2 J (2vr)3 |k-k'|l Ik' 



x 



■ yy ^ + i (k _ k)) . k , r(i)j |k _ k(|) r(i;> |k1) ^ 



(k-k'), (k')j T(r,, |k-k'|) T( fl , |k'| 



We decompose the spatial-spatial part of 6y(k) as follows: 



(2.11) 



e|(k) = 6a e L (k) + ( ^ X utJ ) ©r(k), 



(2.12) 
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so that 



e L (A0 



X 



H inf f d 3 k> a(k-k / )g(k / ) 
2 J (2vr)3 | k _ k /|| | k /|| 

-l dT( V ,\lL-y\)dT(v,\V\) , i fk 



(2.13) 



k')-k' T^lk-k'DT^Ik'l) 



and 



6 T (A;) 



H inf 



1 



d 3 k' afk-k'Vk') 



(2tt) 3 |k-k'|§ Ik'll 



x 



-(k-k') + 



3 (k • k 



^2 



k 2 



k 

2 



./ 2 



T^lk-k'DT^Ik'l 



(2.14) 



3. Scalar Density Perturbations in a Radiation-Dominated Universe 

In the previous section we computed the evolution of the stiff source (i.e., the 
linearized Goldstone modes) and the stress-energy source Q^ u resulting from these 
linearized Golstone modes. In this section we compute the scalar density perturba- 
tions generated by this source, relegating many of the computational details and a 
precise description of the notational conventions to Appendix A. 

For computational simplicity, to make the problem as scale free as possible, we 
assume a radiation-dominated universe, with a{rf) = r\ and a single-component fluid 
with = | and work in momentum space. It follows using the Green functions 
derived in Appendix A that for the Newtonian potential 



(3.1) 



where 



^(M, j) = -(47TC7) [ J df j ff y 1 (kf l /VS) 

o 



x 



e L (k, ff)- 1 - e 00 (k, ff) + ^e T (k, ff) - 



(3.2) 
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and similarly 



i/>(k,ri,y) = +{4irG) (J^ j dfj ff ji(kfj/VS) 

(3.3) 



x 



e L (k, ff) - ie 00 (k, fj) + ^e r (k, 77) - 2e ^ ,j?) 



We now formulate in more quantitative terms the issue of phase dispersion dis- 
cussed in the introduction. We compute the covariance matrix 



(ip(k,r],j) il>(-k,r),j)) (il>(k,r),j) V(-M,y)) 
Cov(k,7/) = ( ) . (3.4) 

(^(k, 77, y) V(-k, 77, j)> (^(k, r/, 1/) ^(-k, 77, y)) 



Here j and y refer to the two types of Bessel functions — or equivalently, the growing 
and decaying modes, respectively. 

The eigenvalues of the covariance matrix provide a measure of the phase disper- 
sion. For the usual adiabatic perturbations from inflation (for which only a growing 

(I ()\ 

mode is present), the covariance matrix is proportional to I ^ I . In the most 

general case with no phase dispersion, this matrix has one vanishing eigenvalue, in- 
dicating a definite fixed linear relationship between ip(k,j) and ip(k,y). By contrast, 

.10' 

a completely random phase gives a covariance matrix proportional to 
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We now compute the matrix element 
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(ijj(k,r],y) ijj(-k',r]',y)) 

,kk' f d s h f d s h 



= Ht nf (4nG)^ ■ J ^ J±]* Jdr ]1 r ] f J d m 4 h(k m /VS) h{k m /V?>) 



(a(k - ki) a(ki) a(-k' + k 2 ) a(-k 2 )) 



x 



|k-ki| 3 /2|ki| 3 / 2 |k'-k2| 3 / 2 |k 2 | 3 / 2 

1 gjp(]k--ki|?7i) djo(kir]i) 
3 <9?7i <9?7i 

f(k-ki) 2 (k-ki) 2 ;2 1 . .„ , . . . . 
+ ( k 2 ~ ^-3-^ " 3^1 Wl k - k il^) Jo(Mi) 



X 



1 gjo(|k^ - k 2 |?? 2 ) dj (k 2 r]2) 
3 9r/ 2 (9r/ 2 
,2 



+ £ {^i + <* ■ >*) - §t} ^{* (|k ' " k2l " 2) » ( ^»} 

integrand with 0^(k) -> 0^(k = 0), B^k') -> 0^(k' = 0) j . 

(3.5) 

The matrix element 

(a(k - ki) a(ki) a(-k' + k 2 ) a(-k 2 )> (3.6) 

as a consequence of Wick's theorem may be decomposed into a sum of products of 
two-point functions with three terms, one of which vanishes, specifically the one that 
would result from computing (^(k, 77 , j ) ) (■?/>(— k', 77, j)). Therefore eqn. (3.6) equals 



(27r) 6 - ( 5 3 (k-k / ) 



5 3 (k 1 -k 2 ) + 5 3 (k-k 1 -k 2 ) 



(3.7) 



allowing us to set k = k' and to consider two cases: (1) k 2 = ki, and (2) k 2 = 
k — ki. By a simple symmetry the two cases give the same contribution, so that eqn. 
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(3.5) becomes 

Myy^TH,^ -5 3 (k-k') 

,k 2 



x 



Hf nf {^G) 2k --5\^-V) Jd 3 h Jd m nf J d m 4 

I9j (|k-ki|?7i) dj (kmi) 



3 3i(km/V3) ji(kr] 2 /V3) 



|k-ki| 3 |ki| 3 



+ 



k i f Jo(|k - ki|?7i) jo(kir}i) 



+ 



3 dr]i drji 

(k-ki) 2 _ k_Jci 2 

k 2 ~ 3 

1 dj (\k - ki|?72) djpjhm) 
3 <9r?2 <9?72 

(k-ki) 2 k-ki 2 



jo(|k - ki|?7i) jo(kir)i) 



k 2 
k 2 \2 



3 ^i !> Jo(|k'-ki|?72) Mhm) 



3(k-ki) 2 l 1 9 



+ 4l^ 2 + (k- ki) - ^^-<| io(|k - ki|i») jo(^ir7 2 ) 



— [ integrand with 0^(k) — > 0^,(k = 0) 

(3.8) 

In order to compare the matrix elements in the most meaningful way, we consider 
the quantities 



F jj (k, V ) = k 3 
F jy (k, V ) = k 3 

Fyy(k,r])=k 3 



ji(kr}/V3) 


ji(kr}/V3) 


V 


V 


jiikrj/VS) 


yi(krj/V3) 


V 


V 


yi(kv/Vz) 


yi(kv/V3) 


V 


V 



Mjj(k;r],r]), 
M jy {k;r],r]), 
■ M yy (k;r],r]). 



(3.9) 



In Figure 1 the F matrix elements are plotted for various k at r] — 1. For kr] ^ 1 the 
decaying modes of the Goldstone field become important, and our calculation here 
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becomes unreliable because of the neglect of these modes. [Because of the scaling 
properties of a radiation-dominated universe, any equal time matrix elements may be 
obtained by a trivial rescaling of equal time matrix elements at r\ = 1.] The /c 3 factor 
has been introduced to make F scale free — for a Harrison-Zeldovich spectrum F would 
be independent of k. Instead we oberve a linear dependence with k in the F matrix 
elements, reflecting the fact that the radiation fluid/scalar gravity modes are excited 
by the Goldstone field source as horizon crossing is approached. On superhorizon 
scales, as k — > 0, the radiation fluid modes are unexcited. 

The degree of phase dispersion is determined from ratio of the eigenvalues of 



the matrix F 



jj 



jy 



Fyj Fyy 



For (krj) & 1, where our calculation is reliable, F is 



roughly proportional to 



Because of the dominance of the Fjj ma- 



1.0 0.17 
0.17 0.05 

trix element, the eigenvectors very nearly coincide with the 'growing' and 'decaying' 
modes ji(kr]/V3)/r] and y\(kr]/y/3)/r] : respectively, with the 'growing' eigenvalue 
dominating over the 'decaying' eigenvalue by a factor of approximately 49. Thus the 
amplitude of the other mode is approximately 15% that of the dominant mode, which 
almost coincides with the growing mode with very little mixing. 




. 0001 L- 1 



Figure 1. The plotted curves (from top to bottom) represent the matrix elements Fjj, 
Fjy, and Fyy at T] = 1 as a function of wavenumber k. 
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4. Discussion 



We have presented a class of plausible inflationary models with Goldstone modes 
that generate non-Gaussian density perturbations that are comparatively easy to com- 
pute. It is of great interest to determine whether the primordial density perturbations 
of our universe inferred from observation are Gaussian, as predicted by the simplest 
inflationary models, or whether they are of a more complicated, non-Gaussian variety. 
Since Gaussian models comprise an infinitesimal fraction of all possible cosmological 
models, in order to define meaningful tests of Gaussianity to consider such tests in 
the abstract is insufficient. Instead proposed tests of Gaussianity must be judged on 
their ability to rule out plausible non-Gaussian cosmological models. 

In all non-Gaussian models based on symmetry-breaking phase transitions — the 
so-called topological defect models — nonlinearity plays a key role, making the predic- 
tions of such models difficult to compute. By contrast, as we have shown, in many 
cases inflationary Goldstone modes may be linearized with very little error, making 
the computation of iV-point density perturbations a matter of evaluating Feynman 
diagrams. 

One shortcoming of the calculations in this paper, one which we plan to remedy 
in future work, is that we have retained only the Goldstone modes corresponding to 
growing modes on superhorizon scales. Because of the neglect of the decaying modes 
of the underlying Goldstone field, the calculations are reliable only on superhorizon 
scales and cannot be trusted on subhorizon scales. To extend the calculation to 
include the decaying modes should pose no real difficulty, and with this extension it 
should be possible to compute all sorts of correlation functions by evaluating Feynman 
diagrams. One should be able to compute the CMB moments in this way and higher- 
order correlation functions, such as the three-point function and beyond. Another 
future direction involves making sky maps of the CMB anisotropy, a problem for 
which computing N-point functions would not be practical and for which some sort 
of nonperturbative renormalization in position space would be required. 

Within the limitations due to neglecting the decaying modes, we were able to 
determine that, somewhat surprisingly, there is very little phase dispersion. The 
absence of significant phase dispersion suggests that the CMB moments from the 
inflationary linearized Goldstone mode theory should exhibit a series of well-defined 
Doppler peaks. 
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After this work was completed, we became aware of an interesting related pa- 
per' 201 that discusses non-Gaussian perturbations from inflation generated by a mas- 
sive scalar field. 
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APPENDIX A— Green's Functions for a Stiff Source 

We work in Newtonian gauge (which is equivalent to the gauge invariant for- 
malism) 1191 and for the stress-energy in addition to the stiff source described by the 
divergenceless tensor assume a single-component perfect fluid with sound speed 
c s . The line element is ds 2 = a 2 (i])-[r] f j iV + h^ l ,]dx fl dx u and rj^ = diag[— 1, +1, +1, +1]. 

For scalar perturbations, hoo = —20, hoi = 0, and = —2ip 5ij, and the 
linearized scalar Einstein equations read as follows: 



-y • [v 2 v> - + n<j>)\ 



5G) 



5GUT) 



(MG) 



Pc5 + Q 



oh 



(87TC7) ■ [p c (l + w)vu + 6 



I J I 



5G{{L) = -= ■ $ + mi) + H<j)+ (2H + W 



(8nG) ■ [ Pc c 2 6 6f + &i(L)}, 
^•[(V,V>'-^V a )(0- 



-V 2 



si 



{Ala) 
(Alb) 

(Ale) 
(Aid) 



where (L) and (T) denote the longitudinal and transverse scalar parts of the spatial- 
spatial tensors, respectively. Here the dots indicate derivatives with respect to con- 
formal time, Ti = (a/a), p c = (3/87iG)H 2 a~ 2 , and V{ = vu (i.e., v is the potential for 
the scalar part of the velocity field). 
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By adding eqn. (Ala) multiplied by c 2 to eqn. (Ale), one obtains 



i> + (m + 3Hc 2 s )ip + H<p + (2H + n 2 + m 2 c 2 s )<p + -v 2 (</> - v) - c 2 v 2 v> 

= (47rG)-[0 i -c 2 0oo]. 

Here Q{j(L) = 5^9^. Working in momentum space and using eqn. (Aid) to eliminate 
cj) in favor of ip, we may write 



^(k) + 3H(1 + cf) ^(k) + (2H + H 2 + 3H 2 c 2 s ) ^(k) + c^(k) 



(47rG) 



L (k) - c 2 O o(k) 

2 27Y ■ 47i + 27Y 2 (l + 3c 2 ) / 

+ 3 r( k ) " ~gr r(k) ^ - e r (k) 



(A3) 



Here we have adopted the convention 



e,,(T)(k) = | i e T( k). 



(-44) 



For the special case of a completely radiation-dominated universe, so that a (77) = r\ 
and c 2 = ^, eqn. (A3) becomes 



iftk) + -^(k) + ^(k) 
r] 3 



(4ttG) 



e L (k) - ie 00 (k) + ^0 T (k) - ^©r(k) 



(A5) 



which is essentially a Bessel equation of order v = |. The homogeneous ver- 
sion of this equation has the general solution ?/>(k, r]) = r(~ x [ip(k, j) ji(ki]/^/3) + 
■0(k, y) yi(ki] / \^3)~\ . From this homogeneous solution we may construct Green func- 
tions, and eqn. (A5) is solved by 



^(M) =r] 1 [ipCk, r], j) j^kr] / V3) +ip(k,r],y)y 1 (kr]/^)] 



(A6) 
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where j\{x) = sinx/x 2 — cosx/x, y\{x) = —cosx/x 2 — sina;/x, 







x (4ttG) 



L (k, fj) - ie 00 (k, fj) + ^6 T (k, fj) - ^©r(k, fj) 



and 



+k f 
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X (4ttG) 



L (k, fj) - ie o(k, fj) + ^©r(k, fj) ~ ^@r(k, fj) 



20 



